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We present theoretical studies of the nature of the collective plasmon resonances of surfaces upon 
which ordered lattices of spherical metallic particles have been deposited. The collective plasmon 
modes, excited by light incident on the surface, are explored for both square and rectangular lattices 
of particles. The particular resonances excited by an incident beam of light depend on the frequency, 
polarization, and angles of incidence. We show that one can create surfaces for which the polarization 
of the reflected light is frequency dependent. The form of the polarization dependent spectra can 
be tuned by choosing materials and the parameters of the nanoparticle array. 
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I. INTRODUCTION 

It is of great interest to create surfaces whose response 
to incident light can be controlled to produce desirable 
effects of various sorts. This may be accomplished in 
diverse ways. For instance one may modulate the sur- 
face profile. Very striking and remarkable reflectivity 
properties are found in systems containing surfaces whose 
roughness is prepared according to prescriptions provided 
by theoretical simulations pQ. We refer the reader to a 
recent review wherein the scattering of light from rough- 
ened surfaces is discussed in detail [2]- 

In addition to modulating the profile of a surface in a 
controlled fashion, one may also deposit material onto it 
in a structured manner. For instance, to enhance pho- 
tocurrents in solar energy devices, metallic nanoparticles 
are deposited on or near the surface of solar devices [31 H]. 
We remark that the nanoparticle arrays in these systems 
are very highly disordered. On a related note, disor- 
dered structures generated by laser bombardment has 
been shown to modify the color of metallic surfaces in 
a profound way [3] . 

Recently, progress has been made towards using the 
plasmonic properties of nanowires [6l [7] or nanoparti- 
cles [51 |S] as optical polarizers, in a fashion reminiscent 
of wire-grid polarizers. This work promises to make po- 
larizers that are microscopic in size, and which possess 
other interesting properties. For instance, the spectral re- 
sponse of nanoparticles can be tuned according to their 
shape [TO] or internal structure [11] . allowing the fabri- 
cation of optical components with new and interesting 
properties. In all these studies, the properties of struc- 
tures on the nanoscale, notably nanoparticles, lead to 
macroscopic optical effects [T2HH]- 
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In this paper, we explore a very different kind of 
nanoparticle-coated surface, namely, a dielectric sub- 
strate upon which an ordered array of sub-wavelength 
metallic nanoparticles [T3] has been deposited [TBlU8| . 
As we shall see from the calculations presented below, 
such a surface exhibits striking optical properties that 
may be tuned by varying the character of the nanoparti- 
cle array. 

There are two important features of such an array. 
First, let £ be a length that characterizes the size of the 
unit cell, and assume that light of (vacuum) wavelength 
A illuminates the structure. If A/l > 1, a condition satis- 
fied for sub- wavelength arrays, there will only be a single 
reflected specular beam, very much as realized for a per- 
fectly flat surface. For particle arrays with linear dimen- 
sions in the range of the wavelength or larger (I > A), one 
realizes additional reflected beams in the form of Bragg 
waves. As the distance between the nanoparticles de- 
creases into the sub-wavelength regime, the additional 
Bragg waves collapse into evanescent waves confined to 
the near vicinity of the surface, leaving a single specular 
reflected beam, once again very much as for a perfectly 
smooth surface. 

A metallic nanoparticle array also supports plasmonic 
resonances that may be excited by incident light. For 
an isolated nanoparticle of radius a <C A made from a 
"plasmonically active" material such as silver (Ag), the 
optically active plasmonic (Mie) resonance lies in the ul- 
traviolet (at Hui 3.5 eV for Ag, where u> is the angu- 
lar frequency of the incident light). Depositing parti- 
cles onto a substrate will cause inter-nanoparticle and 
particle-substrate interactions that redshifts the Mie res- 
onance so it can lie in the visible, as illustrated by previ- 
ous work; see Fig. 9 in Refs. [19] and [20l 

More important for the present study is the role of in- 
teractions between nanoparticles in a dense array. These 
produce collective plasmonic bands whose dispersion re- 
lations and effects on polarization are controlled by the 
properties of the array. As we shall see, these collec- 
tive plasmon modes can be excited by the incident light, 
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with the consequence that the reflectivity of the surface 
becomes highly dependent on both the frequency, the 
angles of incidence, and the polarization of the incident 
light. If the surface is illuminated with unpolarized light, 
the reflected light will be polarized, but the degree of po- 
larization can be strongly frequency dependent, as illus- 
trated by the calculations presented below. By altering 
the microstructure of the array, one may in principle tune 
the polarization characteristics of the reflected light. The 
purpose of this study is to explore these effects, and to 
demonstrate how such effects can be simulated in a way 
that is not computationally demanding. 

A classic and frequently used model to obtain the opti- 
cal response of nanoparticle patterned surfaces is the po- 
larizable dipole model [21] . In this model, each nanopar- 
ticle is treated as a polarizable point dipole, and the effect 
of the substrate is accounted for via image dipoles. For 
nanospheres with intersphere spacing as small as those 
of our geometries, such an approximation is very poor. 
There are highly localized patches of induced charge 
localized around the points of closest contact between 
neighboring spheres, and also around the closest point of 
contact between the sphere and the substrate. One can 
appreciate this from our previous calculation of the spa- 
tial dependence of enhanced fields near an Ag dimer on a 
dielectric surface, as illustrated in Fig. 10 of Refs. IT9land 
I2TJ1 Also in Fig. 5 of Ref.l2"2"lone sees a detailed presenta- 
tion of the electric field near the point of closest contact 
between two spheres when they are excited at their plas- 
mon resonance. It is essential to include higher order 
multipole moments in the description of interactions be- 
tween nearby nanospheres, and between nanospheres and 
a substrate upon which they are deposited, as we have 
done in the results displayed below. 

A more direct approach is that of the discrete dipole 
approximation (DDA) [23], which can be expanded to 
include retardation and substrate effects [53]. The DDA 
is a direct numerical method, relying on discretizing the 
volume of particles into polarizable dipoles. In contrast, 
the method used in this work relies on multipole expan- 
sion of the quasistatic potential [25] . The advantage of 
this method is that it gives results closer to analytical 
mathematics, as well as being computationally efficient, 
assuming that one does not need to include very high 
multipole orders. 

The organization of this paper is as follows. Sectio n [TT] 
contains the necessary theoretical formalism. In Sec. |III| 
we present our studies of the plasmon collective modes of 
arrays of Ag particles supported by an aluminium oxide 
(alumina, A1 2 3 ) substrate, and in Sec. IV we discuss 
the optical reflectivity of two model systems. Concluding 
remarks are presented in Sec. [V] 



substrate as depicted in Fig. [T] The global Cartesian co- 
ordinate system r = (x, y, z) is chosen such that the plane 
z = coincides with the (flat) surface of the substrate 
that is located in the region z < and characterized by 
the dielectric function e_(w). The ambient (z > 0) is 
assumed to be vacuum, and therefore £+ = 1. 

A set of identical Ag nanospheres of radius a are ar- 
ranged on a regular lattice close to the surface of the 
substrate. For each spherical nanoparticle, we associate 
a position vector R^- = (xij , , h + a) pointing from 
the origin of the (global) coordinate system to the center 
of each particle, where the particles are indexed i,j = 
0, ±1, ±2, . . .. For later convenience, we assume a small 
but finite positive value for the parameter h [Fig. l(b)| . 
The particle i = j = is assumed to be located on the 
z-axis so that Roo = (0,0, a + h). Furthermore, each 
sphere ij has associated with it a (local) coordinate sys- 
tem Sij that has its origin located at the center of that 
sphere and its axes oriented parallel to those of the main 
coordinate system. The position vector in Sij we denote 
by i\ , = (rij-,%,^-)- 

The Ag nanospheres are characterized by the dielec- 
tric function e(w). Corrections to the dielectric function 
of the particles due to, e.g., finite size and temperature 
effects are necessary in order to obtain good agreement 
between theoretical predictions and experimental mea- 
surements [26] . This will not be done here, however, 
since comparison to experimental data will not be our 
main concern. Hence, for reasons of simplicity, we have 
assumed bulk values for all dielectric functions. 

In the following we shall consider both square and rect- 
angular lattice structures. Without loss of generality, the 
coordinate system can be oriented so that the lattice con- 
stants b x and b yi corresponding to the directions x and 
y, respectively, obey b x < b y [see Fig. l(a)|. Here, a caret 



above a vector indicates that it is a unit vector. Further- 
more, since the spheres are non-overlapping, we also have 
that 2a < b x . 

Although in general square lattices are formally a sub- 
set of rectangular lattices, we will here restrict the term 
"rectangular lattice" to lattices for which b x < b y . A 
comparison between square and rectangular lattices will 
provide us with an assessment of the range of electro- 
static coupling between the spheres in the lattice. More- 
over, in the limit b y ^ a, where interactions between the 
spheres in the y direction safely can be neglected, the sys- 
tem essentially consists of non-interacting, parallel, linear 
chains of nanoparticles [27] . Moreover, the word "chains" 
will be used about the lines of spheres parallel to the x 
axis, as b r < b v . 



A. Multipole expansion 



II. THEORY 

The system we study consists of a periodic array of 
non-overlapping Ag nanoparticles supported by an A1 2 3 



In this work we will focus on systems for which b y <C A, 
where the quasistatic approximation applies [25] , Then, 
if all materials are assumed to be non-magnetic, the elec- 
tromagnetic properties of the system is fully described by 
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(a) The nanoparticle lattice seen in the xy plane. 
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Ambient (e+) 



»39 




Substrate (e_) 

(b) A view of the lattice in the xz plane. 

FIG. 1. (a) Nanospheres of radius a are arranged on a 
rectangular lattice with lattice constants b x and b y . The z 
axis points out of the substrate, and we always assume b x < 
by. (b) The ambient-substrate interface is located at z = 0. 
The definitions of d, the effective film thickness; a, the sphere 
radius; b x , the lattice constant along x; and the parameter 
h are indicated. The spheres are all characterized by the 
dielectric function e(w). Note that h is exaggerated for clarity. 

the electrostatic potential, ip(r), satisfying the Laplace 
equation: 

V 2 V>(r) = 0. 

Also, the appropriate boundary conditions [35] on the 
surfaces of each sphere and at the interface between the 
substrate and the ambient have to be fulfilled. By def- 
inition, the electric field can be calculated from E(r) = 
— V?/>( r ) [35]. Finally, we assume that the incident optical 
radiation can be modeled as a spatially uniform electric 
field E of angular frequency lj = 2nc/X. As we are em- 
ploying the quasistatic approximation, u) only appears 
in the frequency dependence of the dielectric functions 
£±{ui) and e(u). 

In what follows, we adapt the multipole expansion for- 
malism presented in detail in Ref. 1191 only modifying it 



where necessary to take into account the symmetries and 
the infinite nature of the lattice. The structure of the lat- 
tice requires the potential to satisfy the Bloch-Floquct 
theorem [55] ■ Let ^ij^ij) denote the scalar potential 
of the particle located at (in the global coordinate 
system) and expressed in terms of its local coordinate 
system, r^. We assume that the potential from each 
nanoparticle, ipij, is identical save for a phase factor, due 
to the phase of the incident electric field: 

^(r ii )=V'oo(r ij -)e ik "- ARi '. (1) 

In writing Eq. ([I]), we have introduced the lattice vector 
AR,j = Rjj — Ron = (ib x ,jb v , 0) describing the periodic 
array of spheres, and kn denotes the component of the 
wave vector k of the incident electric field that is parallel 
to the xy plane: k|| = xk x + yk y . Equation ([I]) signi- 
fies that once the scalar potential of sphere i = j = 
is known, it is essentially known for all spheres of the 
lattice. This is a consequence of the Bloch-Floquet the- 
orem, and the form ([!]) is similar to the tight binding 
description of electron energy bands in solids. The Bloch 
phase factor is also assumed to be present in the poten- 
tial of the corresponding image multipole. In passing, we 
note that Eq. (JlJ also predicts that the potential tpiji^ij) 
is invariant under a replacement of kii by k|| + G mn where 
G mn = 2Tr(m/b x ,n/by,0) denotes a reciprocal lattice 
vector (with m and n integers). This invariance follows 
from the fact that a scalar product of a primitive lattice 
vectors from direct space, and one from reciprocal space, 
equals an integer multiple of 2ir [55J. Hence, for the sake 
of the calculation of the potentials, it suffices to consider 
wave vectors kn in the first Brillouin zone. Moreover, 
since this work considers the limit X^$> b y > b x , one may 
take the limit |k|| | — > when calculating the potentials. 

In the formalism presented in Ref. [THl the total poten- 
tial in the ambient, ip+ (r), is expressed as a superposition 
of the potential corresponding to the background electric 
field (Eo), the potential scattered from each nanoparti- 
cle, and the potential from an image multipole designed 
to take the substrate into account: 



tp + (r) 



r ■ E 



E 

i— — oo j — — oo 
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Here the indices i and j indicate that the quantity tjtjj 
represents contribution to the potential from image mul- 
tipoles associated with the sphere ij, and r is the posi- 
tion vector in the global coordinate system. With the 

shorthand notation J2i m ~ S TO =-i> ^ ne mu ltipole 

expansion we use for the scalar potential is given by |25j : 

^oo(roo) =^2A lm r^~ 1 Y l m (e 00 ,(f)oo), 

I in 
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Here, the symbols Ai m and A, are multipole expansion and its absolute value 



coefficients to be determined. The symbol Yj m refers to 
the spherical harmonic functions as described by Ref. 1251 
By symmetry arguments and by matching the boundary 

conditions at z = [29], the coefficients Ai m are related 

(a) 

to the coefficients A\ m ' by the relation 
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We note that the sum over I for practical reasons has to 
be truncated at a finite value L max , meaning that only 
terms for which / < L max are included in the sum. ft 
should be stressed that local fields may not have con- 
verged fully for this cutoff; however, the observed col- 
lective resonances and reflectivity are less sensitive to 
the details of the local fields. This is especially true 
when considering large arrays of metallic nanoparticles, 
in which the interparticle interactions are stronger than 
the particle-substrate interactions. The potential origi- 
nating from other nanoparticles (and images), ipij, can 
be found from ipoo and Eq. ([!]). In order to determine 
the A[ m coefficients, one forms a linear system of equa- 
tions which couples all multipole orders I. How to form 
and solve this system of equations is discussed in detail 
in Ref. M 

Although we are working with lattices that have dis- 
crete translational symmetry with respect to the lattice 
vector, there is little to be gained by use of a Fourier rep- 
resentation of the relevant lattice sums. In mathematical 
terms, the vector R,^ — Rki pointing from sphere kl to 
image ij does not lie in the xy plane: Its polar angle will 
vary depending on the relative position of the image mul- 
tipoles and the sphere on the z-axis. Hence, the presence 
of the substrate, and thus the image multipoles, leads to 
terms in the lattice sums lacking the symmetry needed 
to use the Fourier transform. For this reason, the nec- 
essary sums over spheres |19j are performed directly in 
real space, up to N = 10 unit cells away from the sphere 
on the z-axis. We stress that these sums include both 
nanoparticles and multipole images out to N unit cells. 

In all of our simulations, we have kept a finite value of 
the parameter h to ensure convergence of the spherical 
harmonics expansion. This is because a singularity arises 
at the bottom of the sphere if we choose h = (30) ■ In the 
simulations we present here, we have elevated the layer 
of spheres off the substrate by the amount h = 0.01a. 
We then find very good convergence in terms of the col- 
lective excitations and reflectivity behavior. It should 
be stressed that global behavior, such as the reflectiv- 
ity of the lattice, does not depend on the details of h, 
even though the convergence of the local fields is more 
sensitive to the parameters L max and h. 

As a proxy for the optical response of nanoparticles or 
clusters of nanoparticles, we use the dipole moment, p, 
or its dimensionless analog defined as (in SI units) [Hi] 



pM 



pM 

a 3 e E ' 



(3a) 



pM = (p f p) 



(3b) 



where the superscript f on a quantity indicates its Her- 
mitian conjugate. The components of the dipole moment 
vector can be found from 1191 
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(4a) 
(4b) 

(4c) 



The ith component of the dipole moment is related to the 
incident electric field (E ) in the vicinity of the particle 
by 



ijEo,j, 



(5) 



where onj denotes the polarizability tensor of one of the 
nanoparticles. 

In passing we note that even though the dipole mo- 
ment only depends on the lowest order (1 = 1) expansion 
coefficients, it carries information on higher order reso- 
nances, as the system of equations for the Ai m coefficients 
couple coefficients of all orders [19] . Moreover, for ease 
of comparison, we will below always refer to the dipole 
moment of a single particle. 



B. Surface reflectivity 

The most readily observable quantity that gives an 
indication of plasmonic activity is the surface reflectiv- 
ity. The reflectivity of the system we consider can be 
calculated via several routes. The first approach that 
will be mentioned is due to Bedeaux and Vlieger [29] 
(see also Refs. [521 I3"TH3"5[) . It introduces an equivalent 
geometry consisting of the same ambient and substrate 
as the original geometry, but without the nanoparticles. 
The influence of the latter is accounted for through effec- 
tive boundary conditions for the electromagnetic field on 
the ambient-substrate interface. These boundary con- 
ditions depend on so-called surface susceptibilities that 
contain, for instance, the effect of the size, shape, as- 
pect ratio, and location of the nanoparticles. The sur- 
face susceptibilities are obtained from the multipole co- 
efficients Ai m (Sec. II A I. The approach of Bedeaux and 



Vlieger (29] [3TJ [33] has proven to produce accurate re- 
sults for the surface reflectivity that compares quanti- 
tatively well to experimental measurements [311 1331 135] . 
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It has been used successfully to interpret and invert ex- 
perimental reflectivity data [26l[35]. Recently, this has 
opened for the possibility of detailed in situ and real time 
studies of the growth of supported nanoparticles during 
deposition |35| . However, the accuracy of the Bedeaux- 
Vlieger method comes at a cost: The approach is some- 
what technical, and may be challenging to grasp for those 
not familiar with it. 

Since our main concern in this study is not quantita- 
tive interpretation of experimental reflectivity measure- 
ments, we will instead follow a simpler and more qualita- 
tive route towards obtaining the surface reflectivity. To 
this end, we define above the substrate a thin film region 
of thickness d [Fig. 1(b)] that contains the nanoparticles, 
and when it is homogenized [36- 38 , an effective medium 
results. The classic Maxwell Garnett theory [36TI40] or 
the Bruggeman theory [41] for such films are essentially 
only sensitive to the volume fraction of particles and the 
(bulk) dielectric functions of the materials of the parti- 
cles and substrate, and not to the shape and environment 
(e.g., the ambient-substrate interface) surrounding the 
particles of the lattice. Instead of employing the Maxwell 
Garnett or Bruggeman theories directly, we construct an 
effective medium theory that does depend on other pa- 
rameters, such as the local environment, via the calcu- 
lated polarizability of the metallic particles. The starting 
assumption of our effective medium model is that the lat- 
tice of nanospheres can be represented as an anisotropic 
thin film of thickness d = 2a + h w 2a [Fig. 1(b)] . 

Due to the symmetry properties of the two-dimensional 
rectangular (or square) lattice, the dielectric tensor ( e ) 
will in our choice of coordinate system be diagonal: 



e- 




(6) 



where e x , e y , and e z are the principal dielectric constants. 

Hence, the polarizability tensor, a, must also be diag- 
onal. The macroscopic polarization, P, defined as the 
total dipole moment per unit volume, is thus given by 



Pi 



3 



b x b v d 4-f 



aijEoj 



(7) 



By inserting Eq. ^ into Eq. ([8]), we obtain 

Di = e Q Eo,i + Pi = £ Eoa 



£o X ( S v + blTd aij ) E °' j "X £ iJ E o,j> 
x y / j—i 



so that we get for the component of the dielectric tensor 

1 



(9) 



The polarizability tensor components, a,-j, that appears 
m Eq. j9j, can be related to the multipole coefficients 
Aim by comparing Eq. ([5| to the multipole expansion of 
ip( r ) ESS- This yields the formulas 



1/2 

1/2 

".«." = -'I 2^ J A W- 

1/2 



(10a) 
(10b) 
(10c) 



and a.jj = for i j. Note that the formula for a xx 
[Eq. (10a)] is only valid for E || x, the formula for a yy 
is only valid for Eo || y, and the formula for a zz is only 
valid for Eo || y. Hence, each component of the polariz- 
ability tensor is determined by one simulation each. This 



is the reason why Eq. ( 10 1 seems to disagree with Eq. Q. 
After determining all components of the dielectric tensor 
describing the effective medium thin film, the reflectivity 
of the system can be readily calculated from standard 
theory [40 . In the following, we will denote the surface 
reflectivity Rp(uj) where the subscript f3 — p, s indicates 
the linear polarization of the incident light. Finally, we 
note that whereas we can assume |k||| =0 for the so- 
lution of the Laplace equation, the lateral wave vector 
enters as 

kii = — sin 6*o (cos sin (^o,0) 
" c 

in the reflectance formulas. This is a consequence of the 
fact that the reflectance in general depends on the angle 
of incidence ED] . 



where Pj denotes the ith component of P, denote the 
polarizability tensor components, and subscript indices 
i = 1,2,3 correspond to subscripts x,y,z, respectively. 
Moreover, b x b y d is the volume of the effective thin film 
covering one unit cell of the 2D lattice. In SI units, the 
displacement field, D, is given by 

D = e E + P. (8) 



III. THE COLLECTIVE PLASMON MODES OF 
NANOPARTICLE ARRAYS 

In order to estimate how the interactions between the 
nanoparticles influence their plasmonic resonances, we 
compare the dimensionless dipole moment [Eq. ([3|] of 
an Ag nanoparticle when it is situated in various en- 
vironments in Fig. [2] The dipole moment of a single 
Ag nanoparticle (with no substrate present) is shown in 
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FIG. 2. The absolute value of the dimensionless dipole moment of a single nanosphere in the following environments: (a) 
"hovering" in vacuum, (b) supported by an A1 2 3 substrate, (c) in an Ag dimer supported by an A1 2 3 substrate and with an 
interparticle center-center distance of 2.2a, (d) in an Ag lattice [Figs. 1(a) and 1(b) with b x = 2.2a and b y = 4b x supported 
by an A1 2 3 substrate. For comparison, a perfect metal sphere (of radius much smaller than A) in vacuum excited by a 
homogeneous Eo has p(u) = 1. The parameters common to all subfigures were h = 0.01a and L max = 30. 



Fig. 2(a) and it exhibits a peak at Huj sa 3.5 eV. This 



is the well-known Mie resonance [42] in the quasistatic 
regime. We note that in obtaining the results of Fig. [2] a 
cutoff of L max = 30 was used in the calculations. More- 
over, the same value of L ma x was assumed in obtaining 
all results that will be presented in this paper. 

When the Ag monomer is supported (h = 0.01a) by 
a semi-infinite A1 2 3 substrate (e_ « 2.76 at 3eV, with 
little dispersion), its dipole moment response is as shown 
in Fig. |2(b)| By comparing Figs. 2(a)| and 2(b) it is 
observed that the presence of the substrate leads to a 
rcdshift of the Mie resonance when the electric field is 
perpendicular to the substrate (E || z), but is almost 
unchanged in the case when the electric field is parallel 
to the substrate (Eo _L z). The Mie resonance splits 
due to the breaking of symmetry caused by the presence 
of a substrate. Note that in this case, full convergence 
was not obtained. Nevertheless, the results are shown for 
easy comparison with the other subfigures of Fig. [2] For 
a detailed discussion of this system, see Ref. [43l 

Next, we consider two supported Ag particles placed in 



a dimer configuration with the distance between the par- 
ticle centers being 2.2a and the dimer axis being oriented 
along the x axis. We also assume that it is "hovering" 
h = 0.01a above an A1 2 3 substrate. We observe that 
the resonance for Eq 1 z splits into two resonances, lo- 
cated at different pho ton e nergy for Eo || x relative to 
the E || y case [Fig. 2(c) . This happens because we 



no longer have rotational symmetry about the z axis. If 
instead the incident electric field is oriented along the 
dimer axis (E || x), one gets a redshift of approximately 

0. 45 eV (relative to a corresponding isolated particle) and 
a significant enhancement of the resonance. The particle- 
substrate interactions are not particularly strong for the 
A1 2 3 substrate we consider; it only leads to a small red- 
shift and broadening of the resonance when the incident 
electric field is directed normal to the surface (Eo || z). 
When the incident electric field is directed normal to 
both the surface of the substrate and the dimer axis, 

1. e., Eo || y, the spectrum looks very similar to the case 
Eo -L z of an isolated nanoparticle supported by the same 
substrate [Fig. |2(b)| . 
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FIG. 3. The absolute value of the dimensionless dipole mo- 
ment for a sphere in a rectangular lattice, with anisotropy 
parameters b y /b x = 1.15 to b y /b x — 4; b x was kept constant 
at 2.2a. For reference, the dipole moment of a sphere in a 
square lattice (b x — b y ) is shown in gray shade. In each plot, 
the exciting electric field is directed along different axes: (a) 
Eo || x, (b) Eo || y, and (c) Eo || z. Note that the scale of 
the second axis in (c) is different from the scales used for the 
equivalent axes of the other two subfigures. 



Finally, we study the case when the neighborhood of 
the particle becomes a two-dimensional rectangular lat- 
tice of (identical) nanoparticles with lattice constants 
b x = 2.2a and b y = Ab x , supported by an A1 2 3 substrate 
(h = 0.01a). In this case, the Mie resonance for Ep || x 
undergoes further redshift and enhancement [Fig. 2(d)] . 
As expected, infinite chains of particles cause stronger 
interactions between the particles, causing more redshift 
and resonance enhancement. When the incident electric 
field is transverse to the chains and to the surface nor- 
mal, i.e., Eo || y, the interactions between the chains of 
particles are rather weak and cause no significant red- 
shift. Finally, for Eo normal to the substrate (Eo || z), 
the particle-substrate interactions mainly serve to kill off 
the resonance. This may be understood by recalling that 



the induced dipole in the particle and the image dipole in 
the substrate will have opposite directions and therefore 
partly cancel the effect of each other. 

By comparing the results for the different ratios b y /b x 
shown in Fig. [2] we see that the local environment can 
significantly alter the optical response of nanoparticles. 
In particular, periodic lattices strongly influence the opti- 
cal response of nanoparticles to an incident electric field. 
In Sec. |IV] we will discuss how this can be observed in a 
reflectivity experiment. 

By studying the dipole moment of a single nanosphere 
in a rectangular lattice, one seeks a better understand- 
ing of how the lattice affects the strength and position of 
particle resonances and under which conditions the rect- 
angular lattice turns into a collection of non-interacting 
chains. Hence, in Fig. [3] we show the (single particle) 
dipole moment as a function of the anisotropy parameter 
b y /b x of the lattice. As this parameter increases, the in- 
tersphere interactions along the y-axis weaken [Fig.|3^b)], 
leading to reduced resonance strength and less redshift 
of the resonance. We observe that from b y /b x = 2.64 
to b y /b x — 4, there is little change in both resonance 
strength and position. Results for intermediate values, 
i.e., 2.64 < b y /b x < 4, are virtually indistinguishable 
from those of b y /b x = 4. From this result, we conclude 
that for b y > 46 x , one may ignore interactions between 
the nanosphere chains altogether (at least for our choice 
of b x ). 

When examining Fig. [3]ja), it is clear that if the elec- 
tric field is directed along x, the change in b y does not 
affect the plasmonic resonances significantly. This is be- 
cause the interactions between neighbouring particles is 
strongest along the direction of the incident field. For 
the case when E || z [Fig.[3]Jc)], however, the anisotropy 
affects the plasmonic resonances by increasing the reso- 
nance strength. This effect is mainly caused by the re- 
duced density of particles on the substrate, as neighbour- 
ing particles cause damping of the resonance. This damp- 
ing effect can be explained in a "hand-waving" fashion 
by considering the fact that side-by-side parallel dipolcs 
counteract polarization of each other. 

Simulations have also been run for lattices where the 
anisotropy was kept at b y /b x — 4, but b x (and thus, b y ) 
was increased. As would be expected, increasing the lat- 
tice constants reduces interparticle interactions, and the 
dipole moment behavior becomes more similar to that of 
Fig. |2(b)| Hence, for lattice constants larger than what 
has been explored here, one can safely assume that a 



lattice with b y /b x > 4 can be approximated as a set of 
non-interacting linear chains of nanoparticles. 

Finally, we have investigated the effects of finite |kj||. 
While the effects of a phase difference between neigh- 
bouring particles can be significant, we do not include 
any results in this paper. The reason is that we are pri- 
marily interested in optical fields, where Ikyl <C ft/by. 
Results of simulations for small, but finite, jky |, indicate 
that the approximation kn = is a good one. 
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IV. REFLECTIVITY OF NANOPARTICLE 
ARRAYS 
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FIG. 4. Reflectivity of a square lattice of spheres sup- 
ported by an A1 2 3 substrate for angles of incidence (a) 
(0 O , <M = (0°,0°) and (b) (0 O ,<M = (45°, 0°). The lattice 



parameters were a = 10 nm and b x 



2.2a. The isotropic 



nature of the square lattice means that the reflectivity is iden- 
tical for p-polarized and s-polarized light at normal incidence. 
The dashed lines shows the reflectivity of a substrate with no 
nanoparticles. 

The reflectivity of an interface is an experimentally 
accessible quantity for probing the nanoparticle sys- 
tem. Optical methods have the advantage of being non- 
destructive, and can be used in situ in various challenging 
environments, e.g., vacuum chambers. For these reasons, 
we have calculated the surface reflectivity of a square and 
a rectangular two-dimensional lattice of Ag nanoparticles 
supported by an A1 2 3 substrate. In these calculations, 
the incident field is assumed to be a plane wave that is 



05 



0.45 


1 1 1 1 


1 1 1 1 


0.40 




■ ■ ■ ■ R, 


0.35 




- - R(°) ~ 


0.30 




- 


0.25 




- 


0.20 




- 


fl K 
U. 10 






0.10 






0.05 






0.00 


1,1, 


i . i 



2.8 3.0 3.2 

hui [eV] 

(a) (0 O , 0o) = (0°,0°) 



3.4 



3.6 



0.45 


1 


1 1 


1 1 1 1 


1 


0.40 






- - 4 0) 




0.35 






■ ■ ■ ■ Rs 




0.30 






- ■ - Rs 0) 




0.25 










0.20 










0.15 










0.10 










0.05 










0.00 


i 


i 


i , i 


i 


2.8 


3.0 

(b) (0o,4 


3.2 3.4 

huj [eV] 

o) = (45°, 0°) 


3.6 



FIG. 5. Reflectivity of a rectangular lattice of spheres 
on top of an A1 2 3 substrate. The lattice parameters were 
a — 10 nm, b x — 2.2a, and b y — 2b x . The dashed lines 
shows the reflectivity of a substrate with no nanoparticles. 
(a) At normal incidence [(8o,</>o) — (0°,0°)] we define p- 
polarized light to be polarized along the ir-axis, and (b) 
(0„, 0o) = (45°, 0°). 



either p- or s-polarized. At normal incidence, we define 
ap-polarized (s-polarized) field by E || x (E || y). The 
lattice constants were b x = b y = 2.2a for the square lat- 
tice and b x — 2.2a, b y = 2b x for the rectangular lattice. 
In both cases, a = 10 nm was the radius of the spheres. 
First we address normal incidence. Figures 



4(a) 



and 5(a) present the calculated reflectivity (Sec. II B) for 
the square and rectangular lattices, respectively, at an- 
gles of incid ence (0 o ,0o) = (0°,0°). The red solid curve 
in Fig. |4(a) depicts the reflectivity from the square lat- 
tice at normal incidence, whereas the dashed line (labeled 
i?(°)) shows the reflectivity from a flat, clean A1 2 3 sur- 
face. From this figure it is observed that there is no dif- 
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FIG. 6. The reflectivity of an isotropic and homogeneous 
Ag film at angles of incidence (a) 9q = 0° and (b) 8o = 45° 
(as the film is isotropic, (j>o is irrelevant). The film thickness 
was assumed to be 4.3 nm, which is equivalent to the mass 
thickness of the rectangular lattice discussed in Fig. [5] 



ference (at normal incidence) between p- and s-polarized 
reflected light, as is to be expected due to the symme- 
try of the lattice. Moreover, the presence of the metallic 
particles cause the reflectivity to increase relative to the 
reflectivity of a clean dielectric substrate. This increase is 
especially pronounced near the redshifted Mie resonance. 

For the rectangular lattice at normal incidence, the re- 
flectivity of p-polarized light looks rather different from 

caused by 

b x ^ by. Since b y = 26 



the reflectivity of s-polarized light [Fig. 5(a) 

was assumed, the particles will 
interact more strongly along the x-direction than along 
the y-direction. Based on what was found in Fig. [3j the 
system can be thought of as a set of weakly interact- 
ing linear chains oriented along the x-axis. Along the 



chain, the interactions are strong (p-polarized light has 
a component along x), whereas between the chains (s- 
polarized light is polarized along y) the interactions are 
much weaker. This can be seen directly in the reflectiv- 
ity curves: For p polarization, a large redshift relative 
to the isolated, single particle Mie resonance is observed 
due (mainly) to interparticle interactions. However, for s 
polarization there is essentially no redshift, meaning that 
the resonance occurs close to the single particle Mie reso- 
nance. Note that the interparticle interactions along the 
chain (x) in p polarization produce a much stronger res- 
onance, and hence higher reflectivity, than the interchain 
interactions in s polarization. 

We now turn to non-normal incidence, (9o,<Po) = 
(45°, 0°), for which results are presented in Figs. |4(b)| 
and |5(b)| In this case, the two polarization states give 
rise to different reflectivity curves even for the square lat- 
tice. From Fig. |4(b)| it is observed that the reflectivity 
for s-polarized light is systematically higher than the re- 
flectivity for p-polarized light. Moreover, the difference 
between the two curves is almost independent of energy, 
and the relative change in the reflectivity as a function 
of energy follows closely that of normal incidence. This 
behavior is due to the difference in reflectivity at non- 
normal incidence for the two linear polarizations, as in- 
dicated by the dashed lines in Figs. |4(b)| and 5(b) 



F or the same angles of incidence [(#o,0o) = (45°, 0°)], 
Fig. 5(b) shows that the reflectivity from the rectangu- 



lar lattice is both quantitatively and qualitatively differ- 
ent from that of the square lattice. As the interactions 
between the particles are now significantly stronger for 
electric fields polarized along the x axis, the p-polarized 
reflectivity still has a significantly redshifted peak. For s- 
polarized light, however, the redshift is negligible. This 
means that this surface acts as a "spectral polarizer" : 
It reflects predominantly p polarized light at around 2.9 
eV, whereas at around 3.5 eV, s-polarized light domi- 
nates. The reflectivity can be increased by simply using 
a greater "mass thickness" (the equivalent thickness of 
a homogeneous Ag thin film) of Ag, i.e., by using larger 
nanoparticles. However, the results will be affected by re- 
tardation effects if the nanoparticle radius increases too 
much. 

For comparison, we present in Fig. [6] the reflectivity 
of a continuous, isotropic, and homogeneous Ag film of 
thickness 4.3 nm. This is the equivalent mass thickness 
of the rectangular lattice with b y = 2b x . By comparing 
the results of Figs. |5(b"J| and[6|it is readily observed that 
the corresponding reflectivity curves are rather different. 
Thus, it is feasible to separate the case of a continu- 
ous thin film from that of isolated metallic island films 
through reflectance measurements. In fact, the reflec- 
tivity of the thin continuous film is much closer to the 
reflectivity of the clean substrate, where few features of 
interest can be observed. 
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V. CONCLUSION 

In this paper, we have investigated the collective 
excitations in square and rectangular two-dimensional 
lattices of Ag nanoparticles supported by a dielectric 
(A1 2 3 ) substrate. In particular, we have established 
that for lattices where the lattice constant is greater than 
approximately 4 sphere diameters (b y /b x > 4), the sys- 
tem can be approximated as a collection of independent, 
non-interacting linear chains. It has also been demon- 
strated that if the incident field is polarized along the 
shortest lattice vector, the anisotropic nature of the lat- 
tice leads to collective resonances that are significantly 
redshifted relative to the single particle Mie resonance. 

We have also presented results showing the reflectiv- 
ity of surfaces patterned by such nanoparticle lattices. If 
the lattice is rectangular, the surfaces exhibit the inter- 
esting property that the reflected light possesses different 
colors in the two linear polarizations (p and s). It is rea- 
sonable to assume that similar behavior will be apparent 
in transmission. For technological applications, the ad- 
vantage of the rectangular lattice configuration is that no 
control is necessary over the nanoparticle orientation (cf. 
anisotropic particles on a square lattice), and that reso- 
nance positions can be tuned via the lattice constants. 

Since rectangular lattices of nanoparticles leave much 



room for tunability of the plasmonic and polarization 
characteristics of an interface, we believe further studies 
are in order to gain more insight. Tunability can further 
be extended by employing other materials; nanoshells, 
which allow for tunability through the core and shell 
radii; or through use of truncated spheres or otherwise 
anisotropic particles. For spherical nanoshells, the the- 
ory sketched in this paper could be applied with relatively 
straightforward modifications. The effects of randomness 
on the properties of such lattices is also a topic of in- 
terest [44] . relevant to experimental conditions in which 
some element of randomness is inevitable. The approach 
here is not very computationally demanding, and should 
thus be suited for exploration of random effects. 
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